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STABILITY  ANALYSIS  OF  A  TENSIONED 
STRING  WITH  PERIODIC  SUPPORTS 


1.  INTRODUCTION 


The  motion  of  tensioned  strings  is  typically  studied  to  understand  the  dynamic  response  of 
more  complicated  systems.  Using  a  number  of  assumptions,  the  underlying  differential  equation 
of  motion  is  typically  a  second-order  wave  equation  written  as  a  function  of  space  and  time.  The 
problem  of  an  unrcinforced  string  is  a  classical  continuous  media  problem  and  has  been 
discussed  by  numerous  authors.1'2'3  The  equation  of  motion  of  a  tensioned  string  with  multiple 
sets  of  supports  has  been  derived.-1  The  tensioned  string  with  periodic  stiffeners  has  been 
analyzed  in  the  frequency  domain  with  a  moving  harmonic  force5  and  a  suddenly  applied 
concentrated  force.6  In  these  papers,  the  stability  of  the  system  is  discussed,  particularly  with 
respect  to  the  value  of  the  stiffeners.  There  is  no  extension  of  the  problem  into  the  wavenumber 
domain.  Other  continuous  systems,  specifically  thin  plates,7'8  thick  plates,6  and  beams111  have 
been  modeled  with  periodic  masses,  dampers,  and/or  springs  attached  to  the  medium.  These 
studies  typically  involve  modeling  a  continuous  system  with  one  or  more  differential  equations, 
using  various  boundary  conditions  to  represent  the  mechanical  elements  attached  to  the  system, 
and  then  deriving  a  solution  technique  to  find  the  corresponding  displacements  of  the  system. 

This  report  presents  a  zero- pole  analysis  of  a  tensioned  string  in  wavenumber-frequency 
space  when  it  is  loaded  with  a  continuous  wavenumber- forcing  function  and  a  discrete  point- 
forcing  function.  The  equations  of  motion  and  corresponding  solution  to  both  mechanical  loads 
are  found  using  a  previously  derived  analytical  method.  These  equations  are  then  transferred 
from  an  infinite  summation  series  into  a  continuous  analytical  expression  using  a  series-to- 
trigonometric  mathematical  formula.  Once  this  is  accomplished,  the  dynamic  response  of  the 
system  is  in  zero-pole  format,  allowing  an  examination  of  the  stability  of  the  system  from  a 
parametric  standpoint.  Additionally,  the  locations  of  the  transfer  function  poles  and  zeros  are 
apparent.  A  numerical  example  is  included  to  illustrate  these  effects.  The  dynamic  response  of 
this  system  is  discussed. 
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2.  EQUATION  OF  MOTION 


The  system  model  is  that  of  an  infinite-length,  tensioned  string  attached  to  periodically 
spaced  discrete  stiffeners,  as  shown  in  figure  I .  The  string  is  under  a  tension  of  T(N),  has  a 
constant  mass  per  unit  length  p(kg/m),  and  an  external  load  per  unit  length  of  f{x,t)  (N/m). 

The  stiffeners  are  equally  spaced  at  a  distance  of  L  (m)  in  the  .^-direction  and  each  has  a  stiffness 
of  K  (N/m).  The  model  uses  the  following  assumptions:  (1)  the  forcing  function  acting  on  the 
string  is  at  a  definite  wavenumber  and  frequency  or  is  a  discrete  point  function  at  a  definite 
frequency,  (2)  motion  is  normal  to  the  string  in  the  transverse  direction  (one-dimensional 
system),  (3)  the  string  has  infinite  spatial  extent  in  the.r-direction,  (4)  the  particle  motion  is 
linear,  (5)  the  string  offers  no  resistance  to  bending,  (6)  the  string  is  perfectly  elastic,  and 
(7)  there  is  no  damping  present  in  the  system. 


Figure  /.  Tensioned,  Reinforced  String  with  Coordinate  System 


The  motion  of  the  system  is  governed  by  the  equation' 


<3"iv(x,/)  d2w(xtt)  „  V?3 
T - : - P - -^5 — - -K  X  w(x,t)d(x  -  nL)  =  j (x,t) , 


dx- 


dt 


«=—  00 


(1) 


where  w(x,t)  is  the  transverse  displacement  of  the  string  (m),  x  is  the  position  on  the  string  (m),  t 
is  time  (s),  and  S(x-  nL)  is  the  spatial  Dirac  delta  function  (1/m), 
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3.  SOLUTION  WITH  A  WAVENUMBER-FORCING  FUNCTION 


The  problem  is  first  solved  with  a  wavenumber-forcing  function,  which  is  a  continuous 
exponential  harmonic  function  in  space  and  time,  written  as 

f{x,t)  -  Fexp(ib:)exp(-ity/) „  (2) 


where  ry  is  frequency  (rad/s),  k  is  wavenumber  with  respect  to  the  .v-axis  (rad/m),  and  F  is  the 
magnitude  of  the  distributed  force  (N/m).  The  solution  to  the  problem  is  written  in  series  form 
where  the  string  displacement  is  equal  to  a  sum  of  unknown  functions  of  wavenumber  and 
frequency  multiplied  by  an  exponential  spatially  harmonic  function  in  the  ^-direction  multiplied 
by  an  exponential  harmonic  function  in  time.  The  displacement  becomes 


w(x,t)  =  ^Wn(k,o))  exp(i£„x)  exp(-i  cot ) , 

n=-a o 


(3) 


where 


k 


n 


(4) 


Using  the  principle  of  virtual  work1 1  results  in  the  governing  equation,  which  is  written  as 


,  2  K *S2?°  \F  m  =  0 

(-W  +  pai-WJk, ..a.)--  n  ’ 

L„,-»  1°  m*° 


K  f  F  m  =  0 

A(km ,  (o)Wm  (km  ,  eu)  -  —  £  Wn  (*„  .^)=n  ’ 


(5) 


(6) 
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where 


A(km,0))  =  -Tk2  +  peo2 , 


(7) 


and  k  is  the  wavenumber  of  excitation  (rad/m).  Solving  for  the  coefficients,  and  noting  that  the 
solution  in  the  wavenumber-frequency  domain  is  the  zero11'  order  coefficient,  yields  the 
displacement  divided  by  the  input  force  as 


jr™ 


-  —  \  Z 

\  /  in  — 


+  A  ’(£,<y). 


(8) 


which,  after  some  rearranging,  becomes 


fVo 


(-Tk2  +  pa>2) 


(K/L) 


«-+oc 


-+1-  z 


(K/L) 


(-Tk  2  +  pio1 )  (-Tk2  +  p(o~ ) 


H  =  +CC 


Z 


(K/L) 


»=-^(-Tkl  +  pco  2 ) 


(9) 
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4.  SOLUTION  WITH  A  POINT-FORCING  FUNCTION 


The  problem  is  next  solved  with  a  point-forcing  function,  which  is  a  continuous  exponential 
harmonic  function  in  time  and  a  delta  function  in  space,  written  as 


fix,!)  =  FS{x-x0)cxp{-\cot) , 


(10) 


where  .v0  is  the  spatial  location  of  the  point  force  (m)  and  Fis  the  magnitude  of  the  point  force 

(N).  To  differentiate  this  solution  from  the  solution  found  in  section  3,  the  displacement  variable 
is  changed  from  to  u(xj).  The  system  equation  in  the  wavenumber-frequency  domain, 
when  ,v()  =  0 ,  becomes 


IS  JJ=+O0  £7 

(~k2T  +  pco2)U(k,co)  =  -  ££/(*„.»)  +  -. 

^  n=-o o  ^ 


(ii) 


This  equation  is  solved  in  the  wavenumber  domain  using  previously  developed  analytical 
methods,7'*  and  the  result  is 


U{k,co) 

(F/L) 


{-Tk2  +  peo2) 


l 


-z 


(K/L) 


(~Tk2  +  pco2 ) 


(12) 
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5.  TRANSFORMATION  INTO  ZERO-POLE  FORM 


To  facilitate  an  understanding  of  this  system,  it  is  desirable  to  transform  equations  (9)  and 
(12)  into  a  single  analytical  expression  that  does  not  contain  summations.  This  transformation  is 
accomplished  by  rewriting  the  infinite  series  solution  term  as': 


where  c  =  -Jf  Tp  . 

Inserting  equation  ( 1 3)  into  equation  (9)  results  in 

WoikM  zw{ktm) 

F  p(k,m)  ’ 

where  zw{k,co)  is  an  expression  that  contains  the  system  zeros  for  a  continuous  forcing  function 
and  p{k,co)  is  an  expression  that  contains  the  system  poles  in  the  wavenumber-frequency  plane. 
These  expressions  are 


zw(k,co)  =  Icoyjfp 


K 


L(-Tk2  +pa>2) 


+ 1 


— 1 

o 

>> 

.  2 

r  (oL\ 

+  K  cos 

«)L N 

i(oL) 

sin  — 

-sin 

— 

sin 

l  2  , 

1  2c  )\ 

.  2c  J 

K  2  C  J 

(15) 


and 


p(k,o>)  =  (-Tk2  +  pro2  )\(2co-Jfp) 


sin 


'  kL 

l  2  ) 


-sin 


fcaL' 
2c  j 


+  K  cos  — 

l  2c 


( (oL\  .  ( (oL 


sin  * 

l  2c 


(16) 


Inserting  equation  (13)  into  equation  (12)  results  in 


9 


(17) 


UjkM  _ 

(F/L)  p(k,co)  ' 


where  zu(k,ct))  is  an  expression  that  contains  the  system  zeros  for  a  point-forcing  function  and 

p{kt(o)  is  an  expression  that  contains  the  system  poles  in  the  wavenumber-frequency  plane. 
The  expression  for  zit(k,ot)  is 


zu(k,eo)  ~  2co^fp 


kL 

_  j 

coL) 

1 _ 

( 

c  i  n 

u 

2c  j 

MU 

\ 

kL 

2 


(oL ' 

2 7; 


(18) 


The  system  pole  expression  for  the  discrete  point-loaded  string  is  identical  to  the  system  pole 
expression  for  the  continuous  wavenumber-loaded  string. 

The  terms  for  the  poles  are  now  examined.  The  denominator  in  equation  ( 1 4)  has  two 
distinct  terms  and  both  are  set  equal  to  zero.  The  first  term  is 

-Tk2+pa>2  =  0,  (19) 


and  this  corresponds  to  the  pole  location  of  the  string  without  stiffeners.  This  is  referred  to  as  the 
fundamental  pole  or  resonance  of  the  (unstiffened)  string.  When  K  is  not  equal  to  zero,  the 
numerator  and  the  denominator  in  equations  ( 14)  and  (17)  both  approach  zero  when  the 
relationship  between  wavenumber  and  frequency  in  equation  ( 19)  is  satisfied  and,  thus,  the  term 
given  by  equation  ( 1 9)  does  not  correspond  to  a  pole  of  the  stiffened  system.  The  second  term  is 

(2COyjTp) 

which  can  be  rewritten  as 


•  2 (kL}  .  2 

sin  —  -sin 

l  2  J 


+  K  cos 


coL 

Tc 


.  (  (oL  \ 
sin  — 
l  2  cj 


=  0, 


(20) 


sin' 


coL 

2c 


K 


2COyjTp 


cos 


(£) 


sin 


2  c) 


(21) 
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Note  that  for  the  equality  in  equation  (21 )  to  hold  true,  it  is  required  that 


0  <  sin2 


(  (oL  ^ 

K  ( 

coL  \  . 

cm 

f  coL '' 

Vk) 

t  FT C0S 
2Ct)yjTp  \ 

-  Mil 

,  2  c  J 

\  2  c  j 

(22) 


and  this  occurs  when 


( eoL  mx' 1  ^  K 

tan - -  < - ==, 

l,  2c  2  )  2<OyJTp 


where  m  is  the  largest  integer,  so  that 


0S--  —  Z- 
2c  2  2 


(23) 


(24) 


As  a  result,  the  inequality  given  by  equation  (23)  is  a  necessary  condition  for  poles  to  exist 
in  the  wavenumber-frequency  plane  for  the  stiffened  system.  This  expression  is  not  a  function  of 
wavenumber.  Figure  2  is  a  plot  of  equation  (23)  where  the  ,r-axis  is  the  term  K  /(2 co^Jfp)  and 
they-axis  is  nondimensional  frequency  c oL/2c .  The  gray  area  of  the  plot  is  a  region  that 
corresponds  to  a  location  where  poles  may  be  present  in  the  system  response  and  the  white  area 
is  the  region  where  poles  cannot  exist.  When  the  inequality  in  equation  (23)  is  satisfied,  the 
poles  of  the  stiffened  system  will  reside  at 


arc  sin 


+  mr 


n  =  3,  -2,  —  1,0, 1,2,3,- 


(25) 


where 


6  =  sin2 


f(oL\ 

1 2c  J 


K  {(oL\  .  i(oL\ 

- ■=  cos  —  sin  —  , 

2(o^jfp  V  2c  J  \2c  ) 


(26) 


It 


Nondimensional  Frequency  to  L  /  2  c 


Constant  K  /  [  2  to  (Tp)(1/2)  ] 


Figure  2.  Regions  of  Pole  Existence  (Gray)  and  Nonexistence  (White) 
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and  ktf  is  the  location  of  the  n  indexed  poles  (rad/m). 


The  zeros  of  the  system  with  a  wave  number- forcing  function  are  found  by  setting  the 
numerator  of  equation  ( !4)  equal  to  zero.  This  results  in 


iw(k,eo)  =  Ico-^Tp 


K 


L{-Tk2  +  pat2) 


+  1 


sin ' 


(kL\  .  2 (aL\ 
—  I -sin  “I  — 


V  2c 


f  (oL  \ 
K  cos  — 
l  2c  J 


Si"(^  1  =  0, 


(27) 


which  is  transcendental  in  wavenumber  and  frequency.  The  zeros  of  the  system  with  a  discrete 
forcing  function  are  found  by  setting  the  numerator  of  equation  (17)  equal  to  zero.  This  results 
in 


2h7T 

~T~ 


n  =  -",-3, -2,-1, 1,2,3,” 


(28) 


where  k“  is  the  location  of  the  n  indexed  zeros  (rad/m)  for  the  system  with  point  forcing.  The 

value  of  77 =  0  is  not  a  zero  because  the  fundamental  pole  exists  at  this  location,  which  negates 
the  effects  of  this  specific  zero.  Note  that  the  zero  location  given  in  equation  (28)  is  independent 
of  the  value  of  the  stiffener  K. 
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6.  NUMERICAL  EXAMPLE 


The  theory  and  equations  that  were  developed  in  the  previous  sections  are  now  illustrated 
using  a  numerical  example.  The  system  parameters  are  as  follows:  the  string  tension  T  is  100  N, 
the  density  per  unit  length  p  is  1  kg/m,  the  spring  constant  K  is  3  x  104  N/m,  the  spring  spacing  L 
is  0.05  m,  and  the  computed  wavespeed  c  is  10  m/s.  Figure  3  is  a  plot  of  the  location  of  the 
system  poles,  plotted  in  black,  in  the  wavenumber-frequency  plane.  This  figure  is  valid  for  both 
the  wavenumber  load  and  the  discrete  load.  The  gray  region  corresponds  to  the  gray  region  from 
figure  2  mapped  into  the  wavenumber-frequency  space.  These  gray  regions  are  where  poles  can 
exist,  and  the  white  areas  are  the  regions  where  poles  cannot  exist.  The  dashed  line  in  the  plot  is 
the  function  k  =  cole  and  corresponds  to  the  pole  location  of  the  unstiffened  (i.e.,  no  supports) 
system.  Figure  4  is  a  plot  of  the  response  of  the  string  versus  wavenumber  and  frequency  with  a 
wavenumber  load.  Figure  5  is  a  constant  frequency  slice  of  figure  4  at  225  FIz  and  corresponds 
to  an  area  where  poles  cannot  exist,  and  figure  6  is  a  constant  frequency  slice  at  275  Hz  and 
corresponds  to  an  area  where  poles  can  exist.  Table  1  lists  the  first  eight  pole  locations  of  the 
system  in  the  positive  wavenumber  space  with  the  value  of  n  used  to  calculate  the  pole  location 
at  a  frequency  of  275  Hz.  Figure  3  and  table  1  were  determined  using  the  positive  values  of 
equation  (24).  The  stability  of  the  system  based  on  the  values  of  the  stiffeners  is  discussed  by 
others,5  although  their  analysis  is  only  in  the  frequency  domain.  Figure  7  is  a  plot  of  the  location 
of  the  system  zeros  for  wavenumber  forcing  in  the  wavenumber-frequency  plane.  These  values 
were  determined  by  applying  a  root  finder  to  equation  (27).  The  roots  that  correspond  to 
k  =  co/c  were  discarded  due  to  the  fundamental  pole  also  residing  at  that  location.  Table  2  lists 
the  first  seven  zero  locations  of  the  system  in  the  positive  wavenumber  space  at  a  frequency  of 
275  Hz.  A  zero  with  a  wavenumber  value  of  127.3  rad/m  exists  at  225  Hz.  Numerical 
simulations  suggest  that  for  wavenumber  forcing,  the  zeros  are  almost  collocated  with  the  poles 
at  all  locations  in  the  wavenumber-frequency  plane  except  around  the  region  where  k  =  ±co/c. 
This  concept  can  be  visualized  by  comparing  figure  3  and  figure  7,  where  the  system  poles  and 
zeros  are  in  almost  the  same  location  except  around  the  region  of  k  =  ±co/c .  Note  that  in  figure 
6  the  poles  and  the  zeros  are  almost  at  the  identical  location  in  wavenumber.  Figure  8  is  a  plot  of 
the  response  of  the  string  versus  wavenumber  and  frequency  with  a  discrete  point  load.  Figure  9 
is  a  constant  frequency  slice  of  figure  8  at  225  Hz  and  figure  10  is  a  constant  frequency  slice  at 
275  Hz.  Figure  1 1  is  a  plot  of  the  location  of  the  system  zeros  for  point  forcing  in  the 
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wavenumber-frequency  plane.  Table  3  lists  of  the  first  seven  zero  locations  of  the  system  in  the 
positive  wavenumber  space  with  the  value  of  n  used  to  calculate  the  pole  location  at  frequencies 
of  225  and  275  Hz.  Figure  1 1  and  table  3  were  determined  using  the  positive  values  of 
equation  (28). 


Wa\enumber  (rad/m) 


Figure  3.  Location  of  the  System  Poles  (Black  Lines)  in  the  Wavenumber-Frequency  Plane 

(Gray  areas  are  regions  where  poles  can  exist.  The  dashed  line  is  the  function  k  =  (o!  c. ) 
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Figure  4.  Response  of  the  String  Versus  Wavenumber  and  Frequency  for  an  Applied 

Wavenumber  Load 

(Units  arc  in  dB  rc  m/(N/m) 
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Magnitude  (dB  re  m/(N/m)J  Magnitude  (dB  re  m/{N/m)) 


Figure  5.  Response  of  the  Wavenumber-Forced  System  at  225  Hz 


Figure  6.  Response  of  the  Wavenumber- Forced  System  at  275  Hz 
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Table  1.  Location  of  Poles  in  Wavenumber  for  a  Frequency  of 275  Hz 


Pole  Location 
(rad/m) 

Value  of  n 

33.3 

0 

92.4 

1 

158.9 

-1 

218.0 

2 

284.6 

-2 

343.7 

3 

410.3 

-3 

469.4 

4 

Frequency  (Hz) 


O1 

0  100  200  300  400  500 

Wa\enumber  (rad/m) 

Figure  7.  Location  of  the  System  Zeros  (Black  Lines)  in  the  Wavenumber-Frequency  Plane 

for  Continuous  Wavenumber  Forcing 

(The  dashed  line  is  the  function  k  =  to! c .) 
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Table  2.  Location  of  Zeros  in  Wavenumber  for  the  System  with  a  Wavenumber  Forcing 


Zero  Location  (rad/m)  for 
275  Hz 


35.8 


90.0 


223.8 


283.6 


345.3 


410.0 


470.4 
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Frequency  (Hz) 


Figure  8.  Response  of  the  String  Versus  Wavenumber  and  Frequency  for  an  Applied  Point  Load 

(Units  are  in  dB  re  m/(N/m) 
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Figure  9.  Response  of  the  Point-Forced  System  at  225  Hz 


Figure  10.  Response  of  the  Point- Forced  System  at  275  Hz 
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Frequency  (Hz) 


500 


Figure  II.  Location  of  the  System  Zeros  (Black  Lines)  in  the  Wavenumber- Frequency  Plane 

for  Discrete  Forcing 

(The  dashed  line  is  the  function  k  -  e oi c.) 
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Table  3.  Location  of  Zeros  in  Wavenumber  for  the  System  with  a  Point  Forcing 


Zero  Location  (rad/m) 
for  225  Hz 

Zero  Location  (rad/m)  or 
275  Hz 

Value  of  n 

15.7 

47.1 

-1 

110.0 

78.5 

-2 

235.6 

204.2 

-3 

267.0 

298.5 

1 

361.3 

329.9 

-4 

392.7 

424.1 

2 

486.9 

455.5 

-5 

7.  CONCLUSIONS 


The  zero-pole  response  of  a  tensioned,  reinforced  string  has  been  derived  in  the 
wavenumber-frequency  domain  for  wavenumber-forcing  and  point-forcing  functions.  The 
stability  of  the  system  for  various  stiffener  values  has  been  demonstrated.  A  numerical  example 
was  included  to  illustrate  where  the  poles  and  zeros  are  located  in  the  wavenumber- frequency 
plane.  It  was  shown  that  there  are  regions  where  the  poles  can  and  cannot  exist.  It  was  also 
shown  that  zeros  for  wavenumber  forcing  are  almost  collocated  with  the  poles,  and  the  zeros  for 
point  forcing  can  exist  almost  everywhere. 
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